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We study the critical behavior of the g-state Potts model with random ferromagnetic couplings. 
Working with the cluster representation the partition sum of the model in the large-g limit is 
dominated by a single graph, the fractal properties of which are related to the critical singularities 
of the random Potts model. The optimization problem of finding the dominant graph, is studied on 
the square lattice by simulated annealing and by a combinatorial algorithm. Critical exponents of 
the magnetization and the correlation length are estimated and conformal predictions are compared 
with numerical results. 
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I. INTRODUCTION 

The effect of quenched disorder at a first-order tran- 
sition point is comparatively less understood than the 
same phenomena at a continuous transition point. In the 
latter case rnLsvance-irrelevance criteria, such as the Har- 
ris criteriontHj can be used to decide upon the stabihty 
of the pure fi:^d-point and also perturbation expansions 
are developedo to treat the effect of weak disorder. If the 
transition in the pure system is of first order, neither a 
general relevance criterion, nor a consistent perturbation 
expansion is known to apply around the discontinuity 
fixed point of the pure model. One remarkable except, 
tion is the stability criterion by Aizeaman and Wehrcl 
(based on aji idea of Imry and WortisO, see also by Hui 
and BerkerQ) , which rigorously states that in two dimen- 
sions (2d) any amount of quenched disorder will soften 
the first-order transition in the pure system into a contin- 
uous one. In 3d the same criterion predicts a cross-over 
phenomenon: generally the transition stays discontinu- 
ous for weak disorder, whereas it tjurns to a second-order 
one for sufficiently strong disorderQ. 

Based on the above rigorous results intensive numer- 
ical work has started to clarify the universality class of 
different disordered models, which have a discontinuous 
transition in their pure form. In 2d most of the work 
has been devoted to the g-state Potts model, for which 
the transition point is known from self-duality also in its 
disordered, versions, and in the pure model exact result 
by Baxtem ensures a first-order transition fcpi? > 4. Al- 
though early MonteXjarlo (MC) simulationslMa left space 
for an interpretatiorilj of a q- independent super- universal 
behavior in raiulom systems, later extensive MCll3 and 
transfer matrixEJ calculations consistently determined q- 
dependent magnetic exponents, whereas the correlation 
length exponent, v, was found to be close to the pure 
Ising value, vj = 1, for all q. 

In the large-q limit thermal fluctuations are reduced 
and as a consequence the pure model is soluble in any 
dimension and a perturbation expansion in powers of 



can be performed. In the same limit for the ran- 
dom model at the phase transition point an effective in- 
terface Hamiltonian has been constructed and mapped 
onto the interface Hamiltonian of the random-field Ising 
modeO. This mapping has then been used to relate 
the phase diagram of the two problems and to deduce 
the tricritical exponents of the RBPM at d > 2 dimen- 
sions. However, in the large-q limit no direct calculation 
to study the critical behavior has yet been performed. 
In 2d the presently known information is obtained via 
extrapolation of the results calculated at finite values 
for q. From these estimates no special type of criti- 
cal behavior is expected in the large-g limit. For ex- 
ample the magnetization scaling dimension, Xm, seems, 
to saturate at a finite, non-trivial limiting valuellj't^ 



lim„ 



i{q) « 0.17 — 0.19. However, at this point one 



should note on the presence of strong (logarithmic) cor- 
rections in the form of 1/lnq, see c.f. Fig. 5 in Ref|lJ. 

In the present paper we are going to perform a direct 
investigation of the critical behavior of the random-bond 
Potts model in the large-q limit. As will be shown, in 
that limit the thermal fluctuations are negligible and the 
calculation of the average thermodynamical and correla- 
tion properties of the model is effectively reduced to an 
optimization problem. Here the competition between or- 
dering effects, originating from a tendency to clustering, 
and disordering effects, due to energy gain from quenched 
disorder, plays an important role in determining the opti- 
mal structure. In two dimensions we perform a numerical 
study based on simulated annealing and a combinatorial 
algorithm, and also conformal aspects of the problem are 
investigated. 

The structure of the paper is the following. In Sect. 2 
we introduce the random cluster representation of the 
Potts model and define the equivalent optimization prob- 
lem emerging in the large-g limit. Results obtained from 
the solution of the optimization problem in different 2d 
geometries are presented in Sec. 3 and discussed in Sec. 4. 



1 



II. CLUSTER REPRESENTATION IN THE 
LARGE-Q LIMIT 

We consider the g-state Potts model on a d-dimcnsional 
hyper-cubic lattice with periodic boundary conditions de- 
fined by the Hamiltonian: 



H 
kf 



^ A'y(5(cr,,crj) , 



(2.1) 



where (Ti are g-state Potts variables {<Ti e {1, . . . , g} lo- 
cated at lattice sites i, the sum goes over all nearest 
neighbor pairs {ij) of the lattice, and Kij > are reduced 
ferromagnetic couplings. The d-dimensional hyper-cubic 
lattice corresponds to a graph G = {V, E), where V is the 
set of vertices, which is identical to the lattice sites, and 
E is the set of edges, which is identical with the bonds 
between neighboring sites on the lattice. In the random 
cluster representationtj the partition sum of the model, 
Z, is expressed as a sum over all subsets U C E oi the 
set of edges (or bonds) as: 



(2.2) 



UCE 



where n{U) denotes the number of connected clusters 
in the subgraph G = (V, U) of G, consisting of all lattice 
sites but the reduced set of bonds in U, and = e^'^ — 1 
is the Mayer-function for the coupling Kij . For the latter 
we use the parameterization: 



(2.3) 



Then the contributions from the different graphs to Z 
are expressed in powers of q: 



UCE 



with 



F{U) = n{U) 



(2.4) 



(2.5) 



In the following we consider the large g-limit {q — *■ oo), 
where the partition sum is dominated by the leading term 
given by the maximum value for F: 



Fo = ma^ucE{FiU)} 



(2.6) 



where —Fq corresponds to the free-energy of the system, 
up to a pre-factor of l/(A:Tlng) = const. Let us de- 
note with Uq the subset of E that gives the optimum 
in (l^ej), i.e. Fq = F{Uo), and with Gq = {V,Ua) the 
corresponding dominant graph. Then the energetic con- 
tribution to —Fq is due to the couplings in the dominant 
graph, whereas the entropic term is related to the num- 
ber of connected parts. In what follows, we use the word 



graph when we mean the subgraph G = iV,U) of G de- 
fined by an edge subset U. 

In the pure system, with Wij = w, the structure of the 
dominant graphs in the different thermodynamic phases 
are trivial. Consider a lattice with N ~ spins with 
fully periodic boundary conditions, the number of bonds 
is dN . Then, in the low-temperature phase with w > Wc 
the fully connected graph (V, E) is the dominant graph, 
thus Fq ^ Fc = [dN{l/d + w) + I]. On the other hand 
in the high-temperature phase, w < Wc, the dominant 
contribution is due to the empty graph (V, 0), with a value 
of Fe = N. At Wc = -1/dN, when F^ = F^, there is 
phase co-existence, which means a sharp phase transition 
even in a finite system in the limit oi q ^ oo. In the 
thermodynamic limit we have Wc ~ 0, and the latent 
heat per site is given by AL/N = 1 in our units. 

Introducing disorder, such that Wij can take randomly 
positive and negative values, the question arises, whether 
this trivial structure of the dominant graph persists at 
the transition point, i.e.: Is there still a co-existence 
between two parts of the graph, one being fully con- 
nected, whereas the other^ia empty? To study this prob- 
lem Cardy and Jacobsenll3 have constructed the inter- 
face Hamiltonian, which is then mapped onto that of the 
RFIM. This has lead to the answer that for c? > 2 the 
effect of small disorder is irrelevant, thus there is still 
phase co-existence and thus the transition is of first or- 
der, whereas in d = 2 the phase co-existence is destroyed 
by any amount of disordjsr, in accordance with Aizenman 
and Wehr exact resultsQ. 

In the following we are going to consider the problem 
in 2d where the dominant graph has a non-trivial struc- 
ture. In particular we study the (fractal) properties of the 
largest connected cluster of Go, denoted by F. In the low- 
temperature phase, T < Tc, F is compact, thus the aver- 
age number of points in F is given by [rtr]av cc N — , 
where L is the linear size of the square lattice and here 
and in the following [. . .]av denotes the average over the 
quenched disorder. In the high-temperature phase, for 
T > Tc, [nr]av stays finite and defines the average corre- 
lation length, ^, through [nr]av '^-^ C^- At the transition 
point the average mass is expected to scale as 



(2.7) 



with a fractal dimension d/ < 2ll3. 

The properties of [nr]av are directly related to the 
asymptotic behavior of the average spin-spin correlation 
function, defined in the large-g limit as 



[G{r)U^ma.,,a,))], 



(2.8) 



where (. . .) denotes the thermal and spatial average over 
all pairs of sites i and j with a distance r. We use the fact 
that correlations between two spins are generally zero, 
unless they belong to the same cluster, when G(r) = 1. 
In the case of T < T^, when F is a spanning cluster 
the probability, Pr{L), that a spin belongs to F is given 
by Pr{L) = [7ir]av/A^, whereas the same probability for 



2 



two spins is Pr{LY . From this follows, that the aver- 
age correlations between two spins separated by a large 
distance of r = L is given by: [C(r)]av — Pr{L)^ = 
([?^r]av/^)^- In the low-temperature phase, T < Tc, 
where the average magnetization, [m]av, is defined as 
)[C(r)]av, we obtain: 



Mav = hm -E^^^i^, T <Tc , 

L^oo L 



(2.9) 



whereas at the critical point the average spin-spin corre- 
lations decay as a power: 



2-df, T = Tc. (2.10) 



Finally, in the high-temperature phase, where the av- 
erage size of r is finite the probability to have a con- 
nected cluster of size r is exponentially small, which 
leads to an average correlation function of the form 
[C(r)]av ~ exp(-r/e), for r > ^ 

In the following we specify the form of the disorder, 
where we make use of the simplification that arises due to 
self-duality that holds under special conditions. Accord- 
ing to the results by Kinzel and DomanyB the random 
model is at the critical point, if the distribution, P{'w), 
of Wij is an even function of w, thus P{w) ~ P{~w). For 
convenience we use the bimodal distribution. 



P{w) = pS{w — w) + (1 — p)S{w + LU) , 



(2.11) 



where the critical point is at p = pc ~ 1/2, whereas the 
reduced temperature, t = [T — Tc)/Tc, can be expressed 



t 



-c^(p-l/2), |t|<l 



(2.12) 



Generally we restrict ourselves to the range of disorder 
parameterized as0<w<l/2. We note that for uj ^ 
one recovers the pure model, whereas for u > 1/2 we 
are in the usual percolation limit. Indeed, for the latter 
range of parameters the dominant graph contains all the 
strong bonds, whereas the weak bonds are all absent. 



III. METHODS AND RESULTS 

According to the results presented in the previous Sec- 
tion the solution of the RBPM in the large-q limit is 
equivalent to an optimizati on p roblem with a non-local 
cost- function given by Eq.(2.£). To find the dominant 
graph of the problem we used standard, approximative 
procedures. Most of the results were obtained by the 
method of simulated annealing, but some calculations 
were performed by an approximative combinatorial opti- 
mization algorithm. 

In the procedure of simulated annealing a hypothetical 
temperature variable, T^, is introduced and, after ther- 
malization, is lowered until the hopefully global minimum 
of the cost-function is reached. In practical applications 



we lowered the temperature as T^, = 1/(t— 0.5), in finite 
time-steps r = 1, 2, ... 60, and checked that the resulting 
configuration does not change after further cooling. At 
a fixed temperature in the thermalization MC steps we 
generally used local rules by creating or removing bonds, 
but sometimes we also considered to move a full line of 
bonds. In order to arrive to the global minimum sev- 
eral different starting configurations are considered (at 
least three, sometimes several hundred), and the best 
final configuration was taken. In the investigations gen- 
erally L X L finite samples with linear size up to L = 24 
were considered and periodic boundary conditions were 
used in both directions. For smaller sizes the averaging 
was usually performed over 10000 samples, whereas for 
larger sizes we used several thousands of realizations. 

Alternatively, for small lo (precisely for lo < 0.25) we 
used a combinatorial optimization algorithm that yields 
a configuration that is close to the optimum but not nec- 
essarily equal to it. Actually the worst case bound for 
the ratio of the value Fq of the optimal solution Uq is 
to the value F{U) configuration U that is found by the 
algorithm is only 2/3, which would be too bad for our 
purposes. However, in typical cases the configurations 
produced by the algorithm are much closer, as we checked 
by comparison with the configurations generated by the 
simulatejd annealing method. The algorithm works as 
foUowsta: 

For all sites i let i^-, ix+, iy- and be its left, 
right, lower and upper neighbor, respectively, and denote 
with [iix-), {iix+), {Hy-) and {iiy+) the bonds (edges) 
between these neighboring sites and i. These constitute 
a minimal set of edges that, when removed from G cut 
the site i from the rest of the graph. Let us denote them 
by 

El := {{iix-), (mx+), (wy-), (My+)} (3-1) 
and their weight 



(3.2) 



The minimum cut between any two pairs of sites, i and j, 
(i.e. the set of edges that has a minimum total weight and 
whose removal from G cuts the graph into two disjoint 
subgraphs, one containing i and one containing j) is then 
given either by Ei or Ej, as long as \wij\ < 1/4, as one 
can easily convince one-selves. 

The idea of the algorithm is as follows: Obviously the 
removal of the edges contained in a minimum cut, like in 
Ei for all i, increases the number of components in the 
graph by one, i.e. one wins one unit in the cost function 
F{U) eq.(2.5). On the other hand one looses w{Ei) units 
and when increasing the number of components of the 
graph G one should keep this weight loss as small as 
possible. Therefore we consider a collection of minimum 
cuts as possible candidates of edge sets to be removed 
from G. 
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Let the edge sets be ordered nondecreasing weight, 
such that w{Ei) < w{E2) < ■ ■ ■ < ^(-£^2) and define 
for aU r = 0, 1, 2, . . . , the edge subsets 



Ur^E\\J E, 



(3.3) 



i.e. Uq = E, and with increasing r successively edge sets 
of non-decreasing weight are substracted from E. When 
doing this initially (i.e. for small r) most of the time a site 
will be isolated that before has been connected to a larger 
cluster and therefore frequently (depending on the weight 
of the substraced edges) F{U) will increase, as desired. 
These are the trial configurations for our optimization 
problem and we take the best solution among them, i.e. 
U* such that MU*) = max{F([/^)|r = 0, 1, . . . , L^}. It 
can be shownN that F{U*)/F{uZ^ > 2/3, where Uop 
is the exact optimal solution of ( |2.6| ). With the combi- 
natorial optimization method we could treat larger finite 
systems (up to 128 x 128), than by simulated anneahng 
and the number of configurations we used were between 
10000 and 1000 for smaller and larger systems, respec- 
tively. 



A. Results at the critical point 

First, we tested the relative accuracy of the two 
methods by comparing the value of the obtained cost- 
functions, Fq, for different finite sizes. As a general ten- 
dency simulated annealing has given lower, thus better 
estimates, but the relative difference for L < 16 was very 
small, less then 0.4%. For the largest system we studied 
by simulated annealing, L — 24, the relative difference 
has increased to about 0.6%. We shall later analyze con- 
sequences of the inaccuracy of the min-cut method in 
the magnetic properties of the RBPM. In the following 
illustration we present results which are obtained by the 
more accurate simulated annealing method. 

Typical optimal configurations for different values of uj 
calculated with the same disorder realization for Wij{= 
±u) are presented in Fig.l. The position of the strong 
bonds (wij — +uj) can be obtained from the optimal 
configuration for u > 1/2, since in percolation only these 
bonds are occupied. As seen in the figure for smaller dis- 
order parameter the optimal graph looks to be more com- 
pact, whereas for stronger lo the optimal configurations 
are very close to each other. This fact is a consequence 
of the presence of a finite length-scale in the problem. As 
shown in the Appendix for small lo the system behaves 
uniformly up to a length-scale, Ic, which is estimated as: 



1 

2uj 



(3.4) 



To observe the true asymptotic behavior in the RBPM 
calculation the system size should be larger than this 
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FIG. 1. Typical optimal configurations for different val- 
ues of u) calculated with the same disorder realization for 

Wij{= ±U}). 

value, L > Ic{lo), therefore we restricted ourselves to not 
too small u values. 

Next we analyze the distribution of the largest con- 
nected cluster, r. Inspecting the structure of a typical 
optimal graph in Fig.l we arrive to the conclusion that 
F is a fractal, so that we take the scaling com bin ation 
nr/L''-^, which corresponds to the form in Eq.( |2.7D . In 
Fig. 2 we present a scaling plot of the reduced cluster 
size distribution, where a data collapse can be obtained 
with a fractal dimension of df « 1.8. 

We note that the points, corresponding to the smallest 
system, deviate more from the hypothetical scaling curve, 
which can be attributed to the effect of the finite length- 
scale, Ic- In the inset of Fig. 2 a similar scaling plot is 
presented in the percolation region, i.e. jJpr uj > 1/2, 
where the fractal dimension of percolationE], dp ~ 91/48 
is used. The scaling curves for w < 1/2 and > 1/2 look 
different: for the RBPM the distribution is broad and 
there is a considerable weight for small clusters, whereas 
for percolation the distribution is single peaked without 
a relevant small cluster contribution. 

Next we calculate the average density of the largest 
connected cluster, [nr]av/L^, from the size dependence of 
which the fractal dimension, d/ in Eq.(2.7) and the mag- 
netization exponent, Xm in Eq.(2.10) follows. In Fig. 3 
we have plotted [nr]av/i^ for different finite sizes in a 
log-log scale, using different values of the disorder pa- 
rameter, Lu. In this figure, besides the results obtained 
by simulated annealing, also points calculated by the 
approximate (min-cut) optimization algorithm are pre- 
sented. As seen the min-cut algorithm works satisfac- 
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FIG. 2. Scaling plot of the reduced size distribution of the 
largest cluster at the critical point of the RBPM at to = 0.4 
for different finite systems. A data collapse is obtained with a 
fractal dimension of df ~ 1.8. In the inset the same quantity 
is plotted for percolation, when uj > 1/2 and dp = 91/48. 







0.2 


0.185(30) 


0.25 


0.188(16) 


0.31 


0.165(15) 


0.4 


0.178(13) 


> 0.5 


0.103(2) 



TABLE I. Scaling exponent, Xm, of the average magneti- 
zation for different disorder parameter, uj. The last row with 
LO > 0.5 corresponds to normal percolation where the exact 
value is = 5/48 = 0.104. 



toiy for small systems, L < 16, when the difference in 
the cost-functions calculated by the two methods is also 
very small. For larger sizes, however, which are beyond 
the possibilities of simulated annealing, the error of the 
optimization algorithm increases. Based on the results 
presented in Fig. 3 the min-cut method tends to generate 
a compact cluster in the large system limit. Therefore we 
used the min-cut method only for limited sizes, which are 
anyhow manageable by the simulated annealing method, 
although with much longer computational time. 

Returning to the average density in Fig. 3 one can ob- 
serve that for the disorder parameter in the RBPM range, 

1. e. < w < 1/2, the points fall on nearly parallel 
straight lines having a slope of —2 + df ~ —0.2, where 
df ~ 1.8 corresponds to the value we used in the scal- 
ing plot of the reduced cluster-size distribution in Fig. 

2. The slope of the same line calculated in the percola- 
tion regime, with to > 0.5 is significantly different, it is 
—2+dp ~ —0.1, where dp is close to the fractal dimension 
of 2d percolation. 

The estimates of the magnetization scaling dimension, 
Xm, in Eg. ( p. 10 ) at different disorder parameter, w, are 
summarized in Table I. 

As seen in Table I the average magnetization exponent. 



FIG. 3. Size dependence of the average density of the 
largest connected cluster at different values of the disorder 
parameter, ui, calculated by simulated annealing and by the 
approximative optimization (min-cut) algorithm. Note that 
the min-cut method has a systematic error for larger sys- 
tems. The slope of the curves, si, for different < u < 0.5 
is approximately identical and indicated by a straight line 
with si = —.2, but this slope differs from that of percolation, 
which corresponds to > 0.5, and the related straight line 
has si — —5/48. Typical error of the simulated annealing 
method is indicated by the error bar, whereas the error for 
percolation is smaller than the size of the symbol. 



Xrm is approximately independent of the disorder param- 
eter for < w < 1/2, and its value is within the range of 
Xm ~ 0.17—0.19. This is in agreement with the estimates 
obtained by extrapolating the results calculated at finite 
(/-sljIIj, thus the two limits seem to be interchangeable. 
The apparent variation of Xm with co can be attributed 
to cross-over effects: at a; = the pure system transition, 
whereas at a; = 1/2 the percolation fixed-point is going 
to perturb the value of effective, finite-size dependent ex- 
ponents. 

The average magnetization exponent, Xm, has been 
calculated by another method, which is based on con- 
formal invariancaij. Here we use the result, that in a 
long strip of width, L^,, and with periodic boundary con- 
ditions the average correlation function along the strip 
decay exponentially: 



[((TjCri+„)]av ~ exp(-u/^L„,) 

where the correlation length, , 
asymptotically behaves as: 



27rXr] 



(3.5) 

for large widths 
(3.6) 



In practical calculations we used strips of widths, L.^j = 
2,3,4 and 5, and with such a lengths that in the cal- 
culated correlation function the exponential decay in 
Eq. ( |3.5[ ) seemed not to change by further increase of the 
length. Generally we went at least up to a length of 64 
sites, which has then limited the available widths, L,,,- 
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— 2 


Lyj 3 


= 4 


Lyj — 5 


0.400 
0.423 
0.452 


0.263(9) 
0.267(1) 
0.266(1) 


0.166(4) 
0.168(5) 
0.170(4) 


0.165(5) 
0.167(2) 
0.169(2) 


0.163(6) 
0.163(6) 
0.163(6) 



TABLE II. Numerical estimates for the average magne- 
tization exponent, Xm, using the correlation length-exponent 
relation in Eq.(p.q) for different widths, L^j. 



The calculated exponents for some values of the disorder 
parameter are given in Table II. 

As seen in Table II the size dependence of Xm is very 
weak for Lyj > 3 and the extrapolated value of Xm — 0.17 
is practically independent of the form of the disorder. 
This estimate is compatible with the previous one ob- 
tained by finite-size scaling. The fact, that this latter 
result lies close to the lower bound of the finite-size scal- 
ing one is probably due to the confluent singularity of 
the percolation fixed-point, which is quite strong in the 
region of lu-s we used in the calculation on strips. 

We have also calculated the central charge of the con- 
formal anomaly, c, from the finite-size correction to the 
free-energy per width: 



/o(i^) = /o((X))- — -fO(L-^) 



with the result: 



c = 0.74(1) 



0.51(1) 
In 2 



This is compatible with previous estimateEl c ~ 0.5/ In 2, 
which is obtained by finite-g extrapolation. 



B. Results outside the critical point 

We close our paper by an investigation of the average 
magnetization, [to(L, i)]av, in the vicinity of the critical 
point. In the scaling region, defined as L\t\'^ = 0(1), 
where v is the critical exponent of average correlations, 
the average magnetization is expected to behave as: 



(3.9) 



where rh{y) is some scaling function. The calculated 
magnetizations at different finite size and temperature 
then should collapse to the same scaling function, pro- 
vided the correct critical exponents, Xm and v are used. 
In Fig. 4 we show the result of such a scaling plot, where 
we used = 1, as found approximately in finite-g cal- 
culations, whereas for Xm we used our previous estimate 
obtained through finite-size scaling at the critical point. 
The data collapse in Fig. 4 is satisfactory, however to ob- 
tain a precise estimate on v one needs to extend the cal- 
culations for larger systems. 
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FIG. 4. Scahng plot of the finite-size average magneti- 
zation in the vicinity of the critical point, for a disorder pa- 
rameter u) = 0.4. The scaling exponents we used here are 
Xm = 0.177 and u — 1. 



IV. DISCUSSION 

In this paper the critical behavior of the Potts model 
with non-frustrated, random couplings is studied in the 
large q- limit. We have shown how the calculation of the 
(■3 free energy and the correlation functions of the RBPM 

can be mapped onto an optimization problem, which is 
then studied by simulated annealing and by an approx- 
imate combinatorial optimization a lgori thm. Working 
with the bimodal distribution in Eq. ( ^.ll| ) our results are 
(3-8) compatible with the RG-phase diagram drawn in Fig. 5. 
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FIG. 5. Schematic KG phase diagram of the 2d RBPM 
with varying strength of bimodal disorder, lj. For details see 
the text. 

The pure systems fixed-point {PURE), situated at 
w = 0, is unstable against any weak disorder, thus the 
critical behavior in the range of < w < 0.5 is controlled 
by the disordered fixed point (DIS) . Our numerical cal- 
culation indeed indicate a universality with respect of the 
strength of disorder. Increasing the disorder over uj = 0.5 
we reach the region of attraction of the normal percola- 
tion, and the corresponding fixed point (PERC) is lo- 
cated at = 00. Our RG phase diagram is completed 
by introducing a repulsive tricritical fixed point, TR, at 
u! = 0.5, which separates the regions of attraction of the 
two non-trivial fixed points, DIS and PERC. The sin- 
gular properties of the TR can be quite unusual, since 
the corresponding optimal graph is highly degenerate: 
the possible configurations include all which interpolate 
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between that of the RBPM and that of normal percola- 
tion. 

The behavior of the system at the fixed-point DIS, 
which is the subject of the present paper, is strongly 
dominated by disorder effects, whereas thermal fluctua- 
tions seem to be negligible. Similar, disorder dominated 
critical behavior occur in random quaiitLim spin chains, 
where analytical results are availa bleO^El, and also in 2d 
random quantum ferromagnetsu. Whether exact results 
can be obtained also for the 2d RBPM in the large-q limit 
will be seen in future research. 
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APPENDIX: LENGTH-SCALE IN THE SMALL 
DISORDER LIMIT 

Here we estimate the size, ^, of a step, which is situated 
at the top of a straight surface of a connected cluster, see 
Fig. 6. Using the bimodal distribution in Eq.(2.11) the 



FIG. 6. A connected cluster with a step of Z-points on the 
top of a straight surface. 



existence of the step is connected to the condition: 

W^]>l, (Al) 



21-1 , 

Y. 0> 



i=l 



where Wi — iuj with the same probability, or equiva- 
lently: 



21-1 



> 



1 

2w 



(A2) 



with Pi = ±1. For large I the probability distribution 
of the sum in the l.h.s. of Eq.(A2) is Gaussian, with a 



variance of \/2l — 1. Consequently the average size of the 
step, Ic, scales with a small lo as 



Ic 



1 

2uj 
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